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Abstract 



In the first part of this paper, we prove local interior and boundary gradient 
estimates for p-harmonic functions on general Riemannian manifolds. With 
these estimates, following the strategy in recent work of R. Moser, we prove 
an existence theorem for weak solutions to the level set formulation of the 
l/H (inverse mean curvature) flow for hypersurfaces in ambient manifolds 
satisfying a sharp volume growth assumption. In the second part of this paper, 
we consider two parabolic analogues of the p-Laplace equation and prove 
sharp Li-Yau type gradient estimates for positive solutions to these equations 
on manifolds of nonnegative Ricci curvature. For one of these equations, 
we also prove an entropy monotonicity formula generalizing an earlier such 
formula of the second author for the linear heat equation. As an application 
of this formula, we show that a complete Riemannian manifold with non- 
negative Ricci curvature and sharp L p -logarithmic Sobolev inequality must 
be isometric to Euclidean space. 



Recently, in [M], an interesting connection between the p-harmonic functions 
and the l/H hypersurface flow (also called the inverse mean curvature flow) was 
established. Let v be a positive p-harmonic function, i.e., a function satisfying 



The objective in [M] is to obtain a weak solution to the l/H flow in the level-set 
formulation of [HI] on a co-compact subdomain f2 of R". Where it is sumciently 
regular, such a solution u will satisfy 



that is, (1.2) with p = 1. Moser's strategy in [M] is to obtain weak solutions u to (1.3) 
as limits of solutions vS v ' to (1.2) as p \ 1. To obtain these solutions « , he first 
uses explicit barriers to solve (1.1) with appropriate Dirichlet boundary conditions 
and applies the correspondence above. The key ingredient of the convergence, and to 
the success of this strategy, is to obtain a gradient estimate on the solutions u — u^ p ' . 
The estimate is aided by the following observation. If 



1. Introduction 





(1.3) 




/= |Vu| 2 = (p-1) 



V«| 2 
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then, expressed in terms of /, (1.2) has the equivalent form 

(1.4) (| - l).f p/2 ~ 2 (X7f, Vu) + f^-'Au = f'\ 

Using this equation, the gradient estimate in [M] is proven from out of a boundary 
estimate, which in turn is derived from a C° estimate by way of explicit barriers, 
certain integral estimates on \Vv\, and a Harnack inequality for p-harmonic func- 
tions on Euclidean spaces. The construction of the barriers in particular rely on 
the underlying manifold being Euclidean. The regularization procedure and also the 
reduction of the convergence argument to a certain gradient estimate, uniform in p 
as p — > 1, have some precedence, for example, in [EI], in other contexts. 

In the first part of this paper, we derive interior and boundary gradient estimates 
on a general Riemannian manifold (M, g) via the gradient estimate technique of 
[ChY, LY] and use it to establish an existence result for the l/if-fiow on a class 
of complete Riemannian manifolds. A new feature of our derivation of the local 
estimate is a nonlinear Bochner type formula relating the nonlinear operator with its 
linearization. 

We first start with an interior/local estimate for positive p-harmonic functions, 
which may be of independent interest. 

Theorem 1.1. Assume that v is a positive p-harmonic function on the ball 
B(xo,R), and that on the ball B(xo,R) the sectional curvature of(M,g), Km > — K 2 . 
Then for any e > 0, 

(1.5) sup \Vu\ 2 < 1} U, n + l)b *A +C(n,K,p,R,e) 
where 

Cn = (»^ + ^)_ + ((| + l) m ax { p-l,l } ), 
C(n,K,p,R,e) = (j^ K > + ^ + P -_2)(n - 1) I + K R 

5max{p — 1, l}(n — 1) 
+ (1 - e)R 2 ' 

Note that & p , n , c Pj „ and C(n, K,p, R, e) all stay finite as p — > 1. Hence Theorem 
1.1 effectively gives an estimate for the gradient of the solution to l/Jf-fiow. Also, 
if v is defined globally, by taking R — > oo, and e — > 0, Theorem 1.1 implies that for 
any positive p-harmonic function v, u = — (p — 1) log v satisfies 

(1.6) |V-u| 2 < (n- l) 2 K 2 . 

The constant (n — l) 2 is sharp in light of the results of [Y, LW] for the case p = 
2. From the proof it is evident that when p = 2, one can relax the assumption 
Km > — K 2 by Ric > — (n — l)^ 2 . A direct consequence of this is that if (M,g) 
is a complete manifold with nonnegative sectional curvature, then any positive p- 
harmonic function must be a constant. In fact, using the gradient technique of this 
paper, one can prove that any positive p-harmonic function must be a constant if M 
has sectional curvature bounded from below and its Ricci curvature is nonnegative. 
However, this result can be obtained by only assuming that the manifold has the 
so-called volume doubling property and satisfies a Poincare type inequality (see, for 
example, [Ho] as well as Section 4 for more details). Hence it holds in particular 
on any Riemannian manifold with nonnegative Ricci curvature. On the other hand, 
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it is not clear whether one can obtain an estimate on |Vu| 2 such as (1.6) under the 
weaker assumption of a Ricci curvature lower bound. 

With the help of the interior gradient estimate, by constructing suitable (local) 
barrier functions we can establish the following boundary estimate: 

For every e > 0, there exists po = p(e) > 1 such that if u satisfies equation (1.2) 
on Q, for some 1 < p < Po, then we have the estimate 

(1.7) \Vu\<H+ + e 

where H denotes the mean curvature of d£l and H + (x) = max {H(x), 0}. 

A similar boundary estimate was first proved in [HI] for solutions to a differ- 
ent equation approximating (1.3) (called elliptic regularization) . Our method is a 
modification of theirs. 

With the help of the interior and boundary estimates above, and following the 
general scheme of [M], one can prove the existence of a proper solution (please see 
Section 4 for the definition) for a class of Riemannian manifolds which include the 
asymptotically Euclidean manifolds considered in [HI]. The following is a special 
case implied by our general existence theorem. 

Theorem 1.2. Let M be a complete Riemannian manifold such that its sectional 
curvature Km{x) > —k(r(x)) for some nonincreasing function k(t) with J °° tk(t) dt < 
co. Let Q, be an end of M. Assume that for some p > 1, 



[ 



t \ V(po-l) 



v(nnB(o,t)) 



dt < oo 



and 

lim sup > — —. = 0. 

°°2r<t<oo V(!in%()) 

Then (1-3) has a proper weak solution u with lim^oo |Vuj(x) = 0. 

The volume growth conditions in the theorem are optimal for the existence of 
proper solutions. This is shown in Section 4. 

After the uniform gradient estimate, the key of the proof of the existence theo- 
rem is to construct certain bounded p-harmonic functions and obtain effective C°- 
estimates of such p-harmonic functions at infinity which hold up as p — > 1. Since we 
need to ensure that the limit as p — > 1 is a nonconstant function which goes to +oo 
at the infinity of manifold, the estimate is somewhat delicate. Here we rely crucially 
on the earlier work of Holopainen [Ho]. One may expect more general existence re- 
sults, for example, for manifolds with a Laplacian whose L 2 -spectrum has a positive 
lower bound [LW]. However, it seems that further work, with perhaps more delicate 
estimates, may be needed. 

In the second part of the paper, which starts from Section 5, we consider some 
nonlinear parabolic equations motivated by Theorem 1.1. First, we prove sharp 
gradient estimates of Li-Yau type for two nonlinear parabolic equations associated 
with (1.1). Both estimates proved are sharp in the case that M has nonnegative 
Ricci curvature and provide nonlinear generalizations of Li-Yau's estimate for the 
heat equation. One of these estimates was obtained earlier in [EV1, EV2] (see also 
[V]) for the special case M = R n , in the study of the regularity of nonnegative weak 
solutions. We refer the readers to Section 5 and 6 for the more detailed discussions 
on these results. 

A little surprisingly, we also obtain an entropy formula for a class of nonlinear 
parabolic equations, generalizing the earlier formula for the linear heat equation in 
[N]. More precisely, we show the result below. 
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Theorem 1.3. Let (M,g) be a complete Riemannian manifold. For any p > 1, 
let v be a positive solution to the equation 

(i.8) = (p - lr 1 div(j v«r 2 v«) 

satisfying J u p_1 d/z = 1. T/ien 



^-W p (v,t) = -tp 



A I 



f p/2 1 X7 i X7ju - -i-ajj 
tp 



+ f p 2 RijUiUj v p 1 dfi 



A 



with u = -(p- ljlogw, / = |Vu| 2 , ay = 9 „ - fEf = ^ fc ^ 3i ^r fci for 

any 2-tensor T where (A 13 ) is the inverse of (ay). The entropy 

W p {v,t)= [ (t\X7<p\ p + <p-n)v p - 1 dp, 

J M 



is defined with v p 1 = n- tT/"/ 2 . 4 ". 1 ?-. , where p* — — £- , and is assumed 

to be finite 



The nonlinear heat equation (1.8) has been the object of some previous study. See, 
for example, [V] for a nice survey on the subject. When p = 1, the above theorem 
limits to the following form. 

Theorem 1.4. Let (M,g) be a complete Riemannian manifold. Letu be a solution 

to 

(1.9) ^ -div(|V«r 1 Vw) + |VM| =0 

satisfying that J M e~ u dfi = 1. Then 
d 



di m 



(u,t) = -tj ^j/ 1/2 V a V^M - jg a/ 3 +^+f 1 RijU i u J ^je u da 



with f = \X7u\ 2 , 2 < a, 13 < n are with respect to the orthonormal frame ei = 
,e 2 , • • -,e„. The entropy 

Wi(u,t) — / (t|Vw| + w — nlogi — log u„ — n) e _u d/i, 
J M 

where u>„ is the volume of unit ball in R n , is assumed to be finite. 

Note that (1.9) is the parabolic equation associated to (1.2), the level-set formula- 
tion of the l/if-fiow equation; it has been studied in [He]. In view of the importance 
of the entropy formula of Perelman [P] , we expect that the above result will also play 
a role in understanding the analytical properties of this nonlinear parabolic equation. 
In fact, as a simple consequence of the entropy formula, one can conclude that on 
manifolds with nonnegative Ricci curvature, any positive ancient solution to (1.8) 
must be constant. Similar to the case of p = 2, the entropy W(ys, t) is closely related 
to the optimal L p logarithmic Sobolev inequality (cf. [DD, G]) of Euclidean space. 
When M = M™, the optimal L p Sobolev logarithmic Sobolev inequality implies that 
W(y>,t) > for any with J v p ~ d/i = 1. As another application, we prove the 
following statement: 

Let (M,g) be a complete Riemannian manifold with nonnegative Ricci curvature. 
Assume that the L p -logarithmic Sobolev inequality (6.22) holds with the sharp con- 
stant on M for some p > 1 . Then M is isometric to R n . 

This generalizes the p = 2 case which was originally proved in [BCL] (see also 
[N] for a different proof via the entropy formula for the linear heat equation). 
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The entropy formula coupled with Ricci flow is of special interest. This will be 
the subject of a forth-coming paper. 

Finally we study the localized version of the sharp gradient estimates of Li-Yau 
type mentioned previously for manifolds with lower bound on the sectional curvature. 
This is carried out in the last section. 

Acknowledgement. We would like to thank B. Chow for conservations, K Ecker 
and T. Ilmanen for their interests and helpful discussions. 

2. The proof of Theorem 1.1 

Let u, f be as above. Assume that / > over some region of M. As in [M], 
define 

= div (/^"^(W)) -p/ p/2_1 (Vti, W). 

Here 

A = id+(p-2) 

which can be checked easily to be nonnegative definite in general and positive definite 
for p > 1. Note that the operator C is the linearized operator of the nonlinear 
equation (1.2). 

The first is a computational lemma. 

Lemma 2.1. 

(2.i) an = (4 + Rim ) + (f-i) iv/i 2 r /2 - 2 . 

Here the Hessian ofu. Rij ts the Ricci curvature of M. 

Proof. Direct calculation shows that 
£(f) = (|-l)r /2 - 2 |V/| 2 +r /2 - 1 A/ + (p-2)Au<Vu,V/)/ p/2 - 2 
+(p-2)(|-l)r/ 2 - 3 <VK,V/) 2 

+(p - 2) (W^r /2 - 2 + f im r' 2 - 2 - (v«, v/) 2 r /2 - 3 ) 

-p/ p/a - 1 (Vu,V/>. 



Using 



Af = 2u i3 + 2(VAu, Vw) + 2R ij u i u j 



and combining terms we have that 

C(f) = r /2 - 1 (24+2(VA u ,Vn)+2 J R lJ ^n J )+(|-l)r/ 2 - 2 |V/| 2 
+(p - 2)An<Vu, V/)/^ 2 - 2 + (p - 2)(| - 2)r/ 2 - 3 <V«, V/) 2 

+(p-2) (u^f^r^ + fimf/ 2 - 2 ) - P / p/2 - 1 (v«,v/>. 

Taking the gradient of both sides of (1.4) and computing its product with Vw, we 
have that 

(| - 1)(| - 2)f / 2 - 3 <V/, Vn) 2 + (| - l)/ p / 2 - 2 (f ijUiUj + Uij f iUj ) 
+(|-l)/ p/2_2 A^V/,V ll ) + / p/2 - 1 (VA u ,Vn) = |r /2 - 1 (V/,Vn). 
Combining the above two, we prove the claimed identity (2.1). q.e.d. 
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Now let r](x) = 8 where 6(t) is a cut-off function such that 6(t) = 1 for 

< t < | and 8(t) = for t > 1. Furthermore, take the derivatives of 6> to satisfy 

< 10 and 0" > — 10$ > —10. Here r(x) denotes the distance from some fixed 
xq. Let Q = qf, which vanishes outside B(xo, R)- At the maximum point of Q, it is 
easy to see that 

(2.2) VQ = (V77)/+(V/)7? = 

and 

> C(Q). 

On the other hand, at the maximum point, 
C(Q) = V div (f/ 2 - 1 AVf)- VP f p/2 -\Vu,Vf} 

+f p/2 - 1 (A(Vf), Vt,} + dW(f p/2 AV V ) - pf p/2 (Vu, Vr?) 



V C(f) - (f + l)f>" ^^'^ + / P/2 div(AV,) - p/i <Vu, V,). 



If 1< p < 2, 



?? 7? 

For p > 2 case we have that 

(ACVtri.Vt?) |V7?| 2 
77 77 
The next lemma estimates div( J 4V»7). 

Lemma 2.2. Assume t/iaf, on B(xo,R), the sectional curvature of(M,g) satisfies 
Km > —if 2 . TTien, ai i/ie maximum point of Q, 

1 1 + KR _ 10 max{p - 1, 1} 

~~R? ^2 



div(AVr7) > -20(n + p-2)^^- ™^ + ( p _ 2 )(Vn, Vt?) 



(2.3) +(p _ 2 )|<5^! _(£_!)«. 

Proof. Direct computation yields 

div(AV,) = A7 /+ (p-2)^ + (p-2)Au^M 

(p 2) ^ V "' V ^^ V "' V?7 ^ + (p 2) Ul3Uinj 

Now using (2.2), (1.4) and that = 2uijUi, we can eliminate Au and V/ to arrive 
at 

div(AV?7) = Ar; + (p-2)^j^ + (p-2)(V M ,Vr,) 
+ (p _ 2) p(VK,Vr ? ) 2 _ ( p_ 1) |Vr ? | 2 



'2/77 '2 ' 77 • 

We only need to estimate the first two terms, for which we compute 

Using the Hessian comparison theorem [GW], which states that rij < 1+ ^ r gij , and 
the Laplacian comparison theorem, we have that 



Aijnj < (n + p — 2) 



1 + Kr 
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Noting that 9' = if r < § , we have that 
A V + (p-2)^i = AijTjij 

> -20(n+p-2) 1+ ^ R -10max{p-l,l}-^. 
Taken together, these estimates prove the lemma. q.e.d. 

Remark 2.3. The above lemma is the only place we need to assume that the 
sectional curvature of M is bounded from below by —K 2 . We expect that by some 
judicious choice of cut-off function, one may be able to prove Theorem 1.1 only 
assuming that the Ricci curvature is bounded from below. 

To prove the theorem, we first estimate £(f) from below. We only need to estimate 
it over the points where / > for our purpose of estimating / from above. Choose 
a local orthonormal frame {d} near any such given point so that at the given point 
Vu = |Vu|ei. Then /i = 2ujiUj = 2unUi and for j > 2, fj = 2ujiUi, which implies 
that 

(2-4) 2u kl = -jfe. 

Now (1.4) becomes 

2^ u n =/- (~ p -wii- 



K 2 

3>2 



Hence 



— l j>2 j>2 



> Mil +2^U^ + -^—^(^Ujjf 



J>2 3>2 



n 2 „ 2 1 »2 1 iV „n' 

—ul 1+ 2^u% + — -f + — T (|-l) 



3>2 J 

2 ,p 1W 2 2 p ^/lWiMn 
(o ^ l)/i w i r/ w ii H 7(7 



n- l v 2 " n- 1 n- l v 2 ' / 

Using (2.4), we can replace all the second derivatives of u and arrive at 

- n-l /+ ra-lU + 2 i )f + 2^f 

i,j=l J j>2 J 

where 



x n- 1 V4 2 / 2 J - 
Hence by (2.1), (2.2) we have that 



n — 1 n — 1 r/2 



fP-AjvW _ 2(n _ 1)jf f 2 r- 1 + ( | - 1)/"- 1 ^. 
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Now combining the previous estimates, we have that 

- ( 2 (n - 1) * 2 + 20(n + p - 3)I±™ + 10™xfr-l,l} ) gP _r 

Since 

Q p ~ 2 ( 2(p / 2 _" 1)a + fe^g) <V«, V,) 2 > - ( 2(P f_" 1 1)2 + Q- 1 ^ 
and 

( ^ - 2) Q^<V«, V,) > - W 

\ n — 1 / n — 1 8e r\ 

with 6 Pi „ = (^—^y- — 2^j, we have that 

> f^rf^CiQ) 

2(l-e) 10 i/ %n(n-l) \ 

( 2 (n - 2 + 20(n + p - 3)i±£5 + 10°>«fa-l,l} ^ gp _r 

where 

'2 (p/2-l) 2 (p-2)p 
n- 1 2 



+ 



((| + l)max{p-l,l}) 



Here we have used < -j^-. Theorem 1.1 then follows easily from the above 

inequality. 

3. Boundary estimate 

Let O C M be an open subset such that fl c is compact and d£l is C°° . Again, for 
p > 1, we consider the p-harmonic equation: 

(3.1) A p (v) := div(|V«| p_2 V«) =0 
Any positive solution v of (3.1) gives rise to a solution u of 

(3.2) A p (it) := div (|V«! P_2 Vw) - |Vit| p = 

via the relationship u = (1 — p) log?;. 

Our primary objective in this section is to prove the following boundary estimate, 
which corresponds to a similar result in [HI]. 

Proposition 3.1. For every e > 0, there exists po = p(e) > 1 such that if u 
satisfies equation (3.2) on f2 for some 1 < p < po, then we have the estimate 

(3.3) |Vw| < H+ + e 

where H denotes the mean curvature of d£l and H+(x) = max {H(x), 0}. 
We begin by recording some simple equations. 
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A p (w) = |Vw| P [Aw + (p - 2) — - — - 

\Vw\ 



Lemma 3.2. Suppose <f> : (0,oo) — > R is a smooth function and r(x) = d(x,xo) 
for some fixed xo G M . Fix p > 1. // £ = cut(:r ), t/ien /or a; G M \ E, we have the 
following formula for w(x) — 4>(r(x)) 

ViVj-w = ^'ViVj-r + (f)"VirVjr 

,if-2 /V.. , n ,(VVw)(Vw,Vw) 

(3.4) = |^'| p ~ 2 (>'Ar + (p- 1)0") 

Now we prove a preliminary estimate that will assist us in the proof of Proposition 
3.1. 

Proposition 3.3. Suppose u is a solution to (3.2) on Q, for some 1 < p < n with 
ii = on dQ.. Define 

P(x) := sup { r > 3y G O c wii/i B(y, r) C Q. c and x G dB(y, r) } , 



and put 



l — min ini(:r) 



i? = min^ inf 8{x),^ 



Choose K > so that Rc > — (n — l)Kg on 

S 2 r ■- { x G M \d(x,n c )<'2R} 

Then 

(3.5) |Vu| < ^jp (l + 2VKR) (l-2^f) _1 

on 9f2. 

We should remark that in the case that M has nonnegative Ricci curvature, the 
above estimate can be sharpened to 

which generalizes the one proved in [M]. 

Proof. Consider the function o : (0, 00) — > E defined by 



a{r) : 
If we put 

f' r 

4>(r) ■— 1 — a / cr(p) dp 
for a > to be determined later, we have 4>(R) = 1 and 

Now let u be a solution of (3.2) for some 1 < p < n. Fix xo G d£l, and choose 
yo G fi c such that B(yo,R) C f2 c and xo G dB(yo,R). Setting r(a) = d(yo,x) and 
w(a;) = (j)(r(x)), we have, from (3.4) and (3.6) 

Ap(u5) = |*T~V (Ar + (p-l)^j 

(3.7) = |0'| P ~V ^Ar- (n- 1) (^/K + i 
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away from the cut locus of yo- Since <j>' < 0, applying the Laplacian comparison 
theorem to (3.7) implies that A p (w) > 0, i.e., that w is a subsolution to (3.1) on 
B(yo, 2R) \ {yo} (recall that 37? < inj(yo))- Now choose a to be 

(3.8) a:=(^j\{p)dp^ 

Then we have w(x) — cj> (27?) = on dB(yo, 27?) and w(x) < 1 on dQ. Extend w to 
a subsolution on all of M \ {yo} by 



w(x) :— 



w{x) if x£ B(y ,27?) \{y } 
if x € B(y ,2R) c . 



Since the function v = exp(u/(l —p)) satisfies v = 1 on dQ and v > in Q, it follows 
from the comparison principle that w < v on Q,. In particular, we have < w < v on 
B(yo, (3/2)7?) nfi, Thus if we define z by z(x) := (1 — p) logw(x), z is a non-negative 
supersolution to (3.2) on the annular region 



A := A R 3R := B (yo, 37?/2) \ B (y , R) , 

satisfying z(:r) > 0. If we extend u to be identically in £l c , we have z > u on cL4, 
and hence in all of A by the comparison principle. In particular, at xo, we have, for 
all smooth paths 7 : (— e, e) — > M with 7(0) = xo and 7'(0) = —u(xo) (where v is the 
unit normal to d£l at xo pointing into fi c ), 

m(7W) - mQeq) _ mQtW) ^ £(7W) _ £(7W) - *(7Qeq)) 

Hence 

\Vu\(xo) < \Vz\(x ). 

It remains, then, just to estimate the gradient of z at x . 
Recall from the definition of z that 

|V,| W _<,-l,iggl 

(3.10) =(p-l)a. 
Integrating by parts, we see 

f 2R a(p) dp = (^-^) (R - 2Ra(2R)) - VK — - pa(p) dp. 



So 



(l + 2 (£zi) ^7?) £ «(p) dP > (l + (^) V^p) .(p)dp 



= « ( ^4 ) (i - M2i?)) • 



Since 
we have 



— p 

n — p 

2cr(27?) < 2~, 



r-'lR 

JR 



(p) dp 



< ^JL^ (1 + 2^ W) (l - 2^f) _1 



(p-l)i? 

and the desired inequality (3.5) follows from (3.10). q.e.d. 
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Now we turn to the proof of Proposition 3.1. We use the device in Lemma 3.4 
of [HI] to obtain a supersolution for the solution u, using the rudimentary control 
near the boundary we have already obtained in the previous Proposition in concert 
with the interior gradient estimates of Theorem 1.1 to ensure that the growth of our 
supersolution is sufficiently rapid. 

Proof of Proposition 3.1. Let H denote the mean curvature of an and set H+ — 
max{_ff, 0}. Choose a smooth function w with w = on an such that 

(3.11) H+<^<H + + e 



dv 



on an. Then |Vw| > and 



Aw (VVw) (Vtu, Vw) 
V^f |Vw| : 
dw 



Ai(u>) = -^-^ - 'J^_ 3 - |Vw| 



< 

on an; indeed, for suitably small 8 > 0, both inequalities hold strictly on Us, which 
we define to be the components of {0 < w < S} containing an. Now introduce 



l-w/5' 

Then also w = on an and w(x) — > oo as x approaches dU& \ dQ.. Moreover, a 
computation shows 

Ai(to) = Ai(w) + ( 1 - \— 2 ) |Vw| , 

V (1 - w/6) J 

so Ai(w) < on Us as well. By Proposition 3.3, and the interior estimates in 
Theorem 1.1, for any precompact open set V containing n c , there exists a constant 
C = C(Km, V) depending on the geometry of an and the lower bound of the sectional 
curvature in V (but independent of p\) such that any u satisfying (3.2) with u = 1 
on an satisfies u(x) < C on V. In particular, if we choose such a C = C(5) for Us, 
we see that w > u for any such u on dUs- More precisely we can define Uc to be the 
components of {0 < w < C} in Us- Clearly i»>iion dUc+i- 
But we may also compute 

A p (w) = \Vw\ p Ai(«>) + (p-l)i - ; 

V |Vw| 

which shows that choosing po sufficiently close to 1, depending on C(5) (and implicitly 

on our original e), as well as bounds on our (fixed) barrier w, we have that A p (w) < 

for all 1 < p < po on Uc+i- Therefore, it follows from the comparison principle for 

equation (3.2) that u < w on Uc+i- Hence, recalling (3.11), we have 

du ^ dw dw ^ 

dv ~ dv dv ~ + 
on an as claimed. q.e.d. 

4. 1/H-flow 

The 1/H (or inverse mean curvature) flow is a parabolic evolution equation for 
hypersurfaces. Given an initial embedding Xq : iV n_1 — > M n , the flow may be 
defined parametrically by Nt = Xt(N) where 

dX t 1 
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Here v denotes the outward normal. Alternatively, the flow may be defined in terms 
of a level set formulation, in which case N t is given by d{u < t} for a u satisfying (1.3). 
A theory of weak solutions to (1.3) was established in [HI], based on a variational 
principle involving the functional 

J u (w;K) = I (|Vw| +w|Vm|) d/i 

J K 

for any precompact subset K C O.; the reader is encouraged to consult this paper for 
further details and the motivation behind this theory. In this section, our interest is 
in the general problem of the existence of weak solutions to (1.3), particularly in the 
case in which iVo is the compact boundary of an end f2 of M. 

Definition 4.1. A function u £ C lo ' c (f2) is called a weak solution of (1.3) if for 
every precompact set K C SI and every w € C^^) with w = u in Q \ K, the 
inequality 

(4.1) Ju(u;K) < J u (w;K) 

holds. A weak solution is called proper if lim^oo u(x) = +oo. 
In [M], the functional 

Jl{w;K) = j (j-\X7w\ p +w\X7u\ p ^ dp, 

was introduced for every precompact set K C O. It was shown that if v is p-harmonic 
and u = — (p — 1) log?;, then 

(4.2) JZ{u;K) < Jl{w;K) 

for every w £ W^j! (f2) satisfying w = u in Q, \ K. By a compactness-type argument 
which originated in [HI], it was also shown (p. 82 of [M]) that once one has a uniform 
estimate on | Vu| that is independent of p as p — > 1, one can obtain a weak solution to 

(1.3) from the limit of u as p — » 1. To ensure that the solution one obtains is proper, 
one needs further estimates. In the case of [M], these are furnished by comparison 
with explicit solutions to (1.2) on R n . In our case, we derive suitable estimates from 
results of [Ho] (see also [LT]). 

First we need to recall the notions of p-Green's functions and p-nonparabolicity. 
Let G{x, z) be the Green's function centered at z G M. One can construct G via a 
compact exhaustion as in [LT] for the case p = 2; this was done in [Ho]. To ensure 
that the Green's function so obtained is positive, one must impose conditions on the 
end Q. For the statement of these conditions, we will need to recall the notion of 
homogeneous ends from [Ho], and the associated volume doubling (VD) condition 
and weak Neumann Poincare (WNP) inequality. 

Let o € M be a fixed point. A manifold M (or an end Q, of M) is said to have 
property (VD) if there exists Ci > such that for any x £ M (or G Q) 

V(x,2r) < CiV(x,r). 

Here V(x,r) denotes the volume of the ball B(x,r). It is said to have property 
(WNP) if there exists C 2 (p) > such that 

L M iu - s| * C2r (nb L (x , 2r) |vu|p ) 

where u — y-A ^ J g , . u. For any subset f2 we define f2(i) = Q \ B(o, t). 



GRADIENT ESTIMATES AND AN ENTROPY FORMULA 



13 



Definition 4.2. Assume that (VD) and (WNP) hold on fi. The end Q is called 
p-homogeneous if there exists C3 such that 

(4.3) V(nnB(o,t))<C 3 V(x,±) 

for every x £ dQ.(t). is called uniformly homogeneous if C*2(p) < C 2 for some C' 2 
independent of p as p — > 1. 

As in [LT], we denote the above volume comparison condition (4.3) by (VC). 

Definition 4.3. A end Q is called p-nonparabolic, for some p > 1, if 

roc / , \ l/(p-l) 

(4 ' 4) A (i^m)) dt< °°- 

Tracing the proof of Proposition 5.7 in [Ho] we have the following estimate on p- 
Green's function. Note that the uniformly homogeneous condition implies a uniform 
Harnack constant for positive p-harmonic functions in the Harnack inequality (2.23) 
of [Ho]. 

Proposition 4.4. Assume that is uniformly homogeneous and p -nonpar abolic. 
Then there exists a positive Green's function G(x, z) on M . Furthermore, there exists 
a positive constant Ca independent of p (as p — > 1) such that 

c roo / . \ l/(p-l) 

for every x € f2(r). 

The factor comes from applying Young's inequality in Lemma 5.6 of [Ho]. 

In [LT] and [Ho], it can be seen that, by solving the Dirichlet problem on a compact 
exhaustion, one can obtain a positive p-harmonic function v with v = 1 on d£l and 
lim^oo v(x) — 0, provided the end is p-nonparabolic. It is easy to see from the proof 
of [Ho] that the above estimate (4.5) also holds for such v. This is what is needed 
for our purposes. Alternatively, one can obtain such an estimate from (4.5), using 
the Green's function as an upper barrier and applying the comparison principle. 

Combining the gradient estimates of the previous sections and the above result 
we have the following existence theorem. 

Theorem 4.5. Let M be a complete Riemannian manifold. Assume that Q is a 
uniformly homogeneous end and is po-nonparabolic for some p > 1. Additionally, 
assume that linv^oo V(r) = where 

= 2r™<oo V(finS(o,t))' 

Then there exists a proper solution u to (1.2) on O, with u = on dQ.. 

Before turning to the proof of the theorem, let us first comment on its assumptions. 
The uniformly homogeneous assumption holds in particular for manifolds of so-called 
asymptotically nonnegative sectional curvature. These are manifolds M for which 
there exists a continuous non-increasing function k(t) : [0, 00) — » [0, 00) such that 
Km > —k(t) and J °° tk(t) dt < 00. The reader is referred to [LT] and [H] for 
the verification of the fact that on such manifolds the first two assumptions of the 
theorem are satisfied. Notice that in this case Theorem 1.1 implies that |Vu|(a-) — » 
as x — > 00. Also notice that if 



V(QnB(o,t)) > St 



l+e 
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for some S > and e > 0, then Q is p-nonparabolic for sufficient small p, and satisfies 
limr^oo V(r) = 0. In particular the above theorem can be applied to the asymptoti- 
cally locally Euclidean (ALE) manifolds, which include the cases considered in [HI]. 
We refer the reader to [LT] and [Ho] for other partially homogenous examples, in- 
cluding the manifolds with finite first Betti number and nonnegative Ricci curvature 
outside a compact set. 

In a sense, the extra volume growth condition (lrnv^oo V(r) = 0) imposed in the 
theorem is optimal. Let (M,g) be Hamilton's 2-dimensional 'cigar' [H], i.e, M = R 2 
and g = ds 2 + tanh 2 sd9 2 . Then M has linear volume growth and nonnegative Ricci 
curvature, and is homogeneous, however 

/ tanhs s 



In 



V tanh f 



is a solution to (f .2) on {u > 0} C M but is not proper. Moreover, we can show that 
in fact M admits no proper solution. For suppose that u is such a solution. Let <p 
be a cut-off function. Integration by parts yields that 



/ (p\Vu\ dn< I \V(p\ d/j,. 

J M J M 



I M J M 

Using the fact that (M, g) is of linear volume growth, by choosing suitable <p we can 
conclude that 



/ 

J M 



f 

JO 



|Vw| dp<C 

I M 

for some C depending on the geometry of M. On the other hand, 

r 2ir 

u(R,6)d6 -> oo 

io 

as R — * oo, since u is proper. In particular, we have that 

/ u 3 (s, 6) d6, ds — » oo. 

But, if we denote the point where s = by o, it is easy to see that there exists C\ 
such that 

u s {s,6)d6 <d I \Vu\dp<dC, 

J B(o,R)\B(o,V) 

which is a contradiction. This shows that the volume condition lim r ^oo V(r) — » 
is necessary for the existence of a proper solution. Note that the example can be 
adapted to any dimension. 

Proof, (of Theorem 4.5.) We first construct p-harmonic functions for p < po 
with v' p ' = I on dQ and v^ p \x) — > as x — * oo. This can be done by solving the 
Dirichlet problem on a compact exhaustion and then taking a limit. In view of the 
regularity result of [L], the p-harmonic function will be C 1 ' . On the set where 
|Vm| 7^ 0, elliptic regularity theory implies that both Theorem 1.1 and the boundary 
estimate of Section 3 are valid and can be applied. The existence of a non-trivial 
limit is then ensured by the estimate below. For p < po it is easy to see that 



a: 



f 

J 1r 



t \ V(P-1) 



v(nnB(o,t)) 

A P0 (r) = J™ ( 



dt} <V(r)*o-i ( Apo (r)y- 



where 

i/(po-i) 

v{nnB(o,t)) ' 
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Hence is also p-nonparabolic for p < po- This ensures that there exist satisfying 

i/(p-i) 

dt. 



viP){x) -(p-iy L x) (v(flnB(M))) 



I2r(x) 

Combining the previous inequalities, we have that for r(x) >> 1 

u (p \x) > 2(p-l)]og(p-l)-(p-l)logC4-(p-l)log(A*(r(x))) 

-?°^£]og(V(r(x))). 

Po - 1 

By the gradient estimates in the previous sections we can conclude that converges 
locally uniformly to a limit function u(x) € C^(Q). Moreover, u(x) satisfies that 

u(x) > - log(V(r(a;))) 

for r(x) >> 1. By the compactness argument from [M], Theorem 1.1, and Proposi- 
tion 3.1, we can conclude that u is a weak proper solution to (1.2). q.e.d. 

5. The p-Laplacian heat equation 

Motivated by Theorem 1.1, we consider smooth solutions to the parabolic analog 
of (1.1), namely 

(5.1) ^ = div(|V«r- 2 Vv) 

for p > 1. This nonlinear evolution equation is the gradient flow for the p-energy 
functional 

E v (v)= f \Vv\ v d[i 

J M 

and has been studied rather extensively - see, for example, [Ba], [V], and the ref- 
erences therein. For a given smooth solution v of (5.1), it will be useful for us to 
consider the linearization of the operator A p at v, given by 

£(V0 = div(ft p/2 " 1 A(V^.)) 
Here h = |Vu| 2 and A is the tensor introduced in Section 2, namely 

The main result in this section is essentially a consequence of the following calcu- 
lation. 

Lemma 5.1. Suppose v : M x [0, T) — > R is a smooth, positive solution to (5.1) 
with p > 1. For any a > 0, define 

F„ 

tr v 

Then, on the region |V«| > 0, the following estimate holds 

2 D . ....... \ ..2 



£-i)F a = P h*>- 2 



dt 



Vij 1 ViVj 



v p — 1 V 2 



A 



(5.2) 

In particular, if He > —Kg for some K > 0, we have 
(5.3) 



^M +(a _ 1)(p _ 2) ^ 



h p/2-l 



+ (p - 2)F X 2 - 2(p - 1)^— — (VF a , Vw). 



-pif^-5--2(p-l)- (VF Q ,Vn). 
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Here, for a two-tensor T, we write \T\\ := A ik A jl T i:j T kl . 

Proof. First we have the following formula resembling Lemma 2.1. 
Lemma 5.2. Let h = \Vv\ 2 . Then 

(5.4) ^ - cjh = -2^ /2 - 1 (|VV^| 2 + RijViVj) - (| - l)/i p/2 " 2 |V/i| 2 . 

The proof of the above identity is a straightforward calculation and very similar 
to that of Lemma 2.1. Hence we leave the details to interested readers. 
By the definition of C, we have that 

(5.5) (£-£)v t = 0, 



9t. 

and with the help of Lemma 5.2, we find that 

(5.6) (E - h p/2 = ph p - 2 [|VVwji + Ravi- 

Using the general formula 

£(—) = —£(Q) - t4tA«) 

\V~t ) V~< 



and the identity 



- - —hf'^VQ, Vv) + 7(7 + l)(p - l)/, p/2 - ■-. 



9 

£ - — ) v = (p - 2)v t 



dt 

together with (5.6) and (5.5), we then obtain 



(5.7) 



and 



(5.8) 



Now, 



~ r) \ h p /' 2 n 



h p ~ 2 h v 
4 P(P - 1) — s-VijViVj + 6(p - 1) — 



h p/2 

+ 2(p-l)^v t . 



V7 T7 u p/2-lVikVi V k VfcWt V t 

V k F a =ph p/ — 5 2^r -a ha^ 

so, after multiplying both sides by a, (5.8) becomes 



+ 4(p-l)^+2(p-l)- (VF a ,V«>. 
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Combining this with (5.7), we have 



d_ 

at 



h p ~ 



\VVv\a + RijViVj] -2(p-2) 



h p ' 2 



-v t 



(5.9) 



h p ~ 2 h p 
- MP - IJ—T-VijViVj + 2(p - l)-j 



+ a (p_2)^-2(p-l) 



ft p/2-l 



Equation (5.2) then follows from (5.9) and the identity 

- 2 va 1 ViVi 2 h p ~ 2 



{VF a ,Vv). 



ph p ~ 



Vi 



v p — 1 V 2 
For (5.3), we observe that 



P- 



Fi=h 



P/2-1 



VVv\ 2 A - 2p(p- V 



2 h p 

— VijViVj+p—. 



(p — 1) V 



ViVj Vij 
V 



where 



so that 



ViVj 



(p-1) V 2 



ph 



p-2 



(p-1) V 2 
2 



ViVj Vij 
V 



(p — 1) V 2 



>P F ? 

n 



by the standard inequality n|T|^ > {%xa T) 2 for a two-tensor T. 



q.e.d. 



In the case a = 1, it is convenient to consider the expression of the above equations 
in terms of the "pressure" quantity 

P — 1 £=2 
= -V P- 1 . 



1 Vt) g Vt; 



(5.10) 
Then 

and 

A p {<j>) = |V0r- 2 tr A (VV^) = -Fi 
so that in the case a = 1, the identity (5.2) has the equivalent form 



__, /VVt> 

VV0 = up- 1 

V v 



(5.11) 



0i 



£ ) A P {<f>) > P \ V0j 2(p - 2) (|VV0& + iii^i^) + (p - 2)A P (0) 2 
-2(p-l)(vA p (0),V^). 



When Rc > 0, (5.3) suggests, by way of the maximum principle, the global esti- 
mate 



K{4>) = ^ 



|Vt;| p _n/3 



(5.12) 

t) v 2 t 

for /3 = l/(p + n(p — 2)) and p > 2n/(n + 1). Such an estimate, analogous to the 
differential Harnack estimate in [LY] for the heat equation, would indeed be sharp 
in view of the explicit source- type solutions (see, e.g., [Ba]) 



H p (x,t) 



1 

t"0 



1 + 



/3p-i (2 -p) f\x\\ " 



p-i 

p-2 



to equation (5.1) on R" for which (5.12) is an equality. However, (5.1) is degenerate 
where |Vv| = 0, and our calculations, carried out in the region |Vi>| > 0, are thus as 
yet insufficient to draw such a conclusion. Following [EV2], we therefore consider a 
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family of strictly parabolic equations which approximate (5.1) and by proving anal- 
ogous estimates for the corresponding quantities of the solutions to the approximate 
equations, we may obtain an estimate of the above form via a limiting procedure. 
The precise statement of the result, which was obtained in [EV2] in the case M = R n 
(albeit for a broader class of solutions) , is the following. We restrict our attention to 
smooth solutions. 

Theorem 5.3. Suppose (M n ,g) is a complete Riemannian manifold with non- 
negative Ricci curvature, and suppose v is a smooth, nonnegative solution to (5.1) 
with p > 2n/(n + 1) for t G [0, fi). Then, for all t € (0, Q), one has 

|V«|" ldv np 

where 

qU if Jfr<P<2 

\(p-l)/3 if P>2 

Remark 5.4. In the case M = R, the above result - with the sharp constant 
(3 = l/2(p - 1) - was obtained in [EV1] for all p > 1. 

Proof. The argument follows almost exactly as in that of Section 2 of [EV2], 
where the equation satisfied by the pressure <j> = (p— l)/(p — 2)v^ p ~ 2 ^^ p ~ 1 \ namely, 

^ = ^^divdv^r^^) + m p , 

is approximated by the equation 

(5.14) ^ = ^-^ E div(^(|V0e|)V0 e ) + Ve(|V0 e |) 

for judicious choices of <p £ and i/> £ depending on e > 0. For completeness, we outline 
the argument below. 

First, take € C°°([0,co)) to satisfy 

Ur) = { p - 1 if ^ where a e = { <<* ~ * P < 2 

SW \l if re[0,a e ] \e(p-l)-V72 if p > 2. 

On the interval (a e ,r), choose £ e to be non-decreasing if p < 2, non- increasing if 
p > 2, and to satisfy |r^(r)| < eCe( r )- Then let <p e satisfy 

r^=C(r)-l, 

so that ip e {r) is constant for small r and equal to r p ~ 2 for r > e. Finally, define 



Mr) ~ " 



—P-j (r 2 <fc{r) - J sip,(s)ds^j . 



With these choices, (5.14) is strictly parabolic and we obtain a smooth solution <j> e 
for each e > 0. 

Next, we derive an equation analogous to (5.11) for Fi, e := — div(y> e (|V0 e |)V0 € ). 
Define 

B ij = a + and B'J = g ij - 

pifie p 2 

where p = |V0 e |. Introducing the operator 

L £ ($) == div (t^l^B.V^ 



GRADIENT ESTIMATES AND AN ENTROPY FORMULA 



19 



we compute 



r e - - J F M = 2<^ e (VFi, e> V0 e ) + J^jCeJi,. 



(5.15) 



+ -riC«^ [|W^j| e +Rc(V0 e , 



+ CM f-^C, + 2^) H 2 + ^«tf(A#, - H), 

\P 1 S,e / P 1 

where 

S = ^S^ViV.-^VfcV^eV^.VVe and = [VV0 e (V<?!) e , V0 e )] 2 . 

Now, F 1>€ = —(pcB^ViVjtjf^, and when p £ [0, o e ] U [e, co), the final five terms on the 
left side of (5.15) reduce to 



jCe^i,. + ^ZTC«V« [|VV0 e || e + Rc(V0 £ ,V<M] 
>-M P -2 + £)f?,. 

p - 1 \ 71/ 



(5-16) > TT^ F M- 

max{l,p — 1 jnp 

When p £ (a E , e), we need to do a little extra estimation. First, note that S > -ff 2 
(so that the fourth term in (5.15) is non- negative for sufficiently small e by the 
assumptions on £' e ). Also, one has 

Fh<—y* (A^-Hf+n^H 2 
n — 1 

(5.17) <n^ 2 |VV0 e || +ng<p 2 e H 2 , 



and similarly for any p > 0, that 

(5.18) H(Acj> e -H) > -^H'- ^— -|VV0. 



f ^_2_i |V v*,|, 



With these inequalities, equation (5.15), and the assumptions on £' e , we can find a 
constant C, independent of e, such that 



^ p4r [(p- 2 ) F ^ + (p- (lW0 e || e + £h 2 )] 



P" 

,2 



(5.19) > ^ + 

max{l,p — 1} \ n J 

Since we have already noted that the stronger inequality (5.16) holds when p < a e 
or p > €, one may appeal to the maximum principle to deduce in general that 

nB 



(l-Ce)t' 



Letting e — > 0, and making use of appropriate energy estimates, as in [EV2], we can 
obtain a limiting 4>, and consequently a solution v to (5.1), for which the inequality 
(5.13) holds. q.e.d. 
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6. Another nonlinear parabolic equation 

Note that with the above parabolic estimate, one can not recover the elliptic result 
proved in Section 2, thus in this section, we consider instead the parabolic equation 
associated to the operator A p . Let u be a solution to the equation 

(6.1) ^ - A p (u) = ^ - div(|V<- 2 Vu) + \Vu\ v = 0. 

This nonlinear parabolic equation has also been studied in the literature. See, for 
example, [V] and the references therein. 

The corresponding equation for v = exp(— ^rr) is 

(6.2) ^^ = (p _i)P-i div (| Vu r 2 V«). 
Recall from Section 2 the operator C defined as 

£(V)=div(r /2 - 1 A(W))- P / p/2 - 1 <Vu,V^), 
which is the linearized operator of A p . Note that 

A-=^ + (p-2 ) ^i=^ + (p-2)^j. 

For any u, letting / = | Vu| 2 , the proof of Lemma 2.1 yields, the nonlinear Bochner 
formula 

(6.3) £(/) = 2f' 2 -\u% + RijUiUj) + 2(X7u, V(A P ( W ))> + (| - l)|V/| 2 r /2 - 2 . 
A bit of computation together with this formula yields the following lemma. 

Lemma 6.1. Let f = |Vitj 2 . Then 

(6.4) (§- t -Au t = 



(6.5) [j t - C )f = -2/" /i, " 1 K-+^Ui«i)-(^-l)/ , ' /2 - 2 |V/| 2 . 
An immediate consequence is the following result. 

Corollary 6.2. For any a > 0, Zei F Q = |Vm| p + au t = f p/2 + au t . Then 



(6.6) 



^ - £^ F a = -p/ p - 2 (|VVu|i + flyuity) ■ 



Notice that (6.1) implies that 

Fi = / p/2 - 1 tr A (VVu). 

Hence 

f-' 2 \Wu\ 2 A > -fI 

n 

In the case that M is compact with nonnegative Ricci curvature, we obtain at once 
the following global estimate. 

Theorem 6.3. Let M be a compact manifold with nonnegative Ricci curvature. 
Let v be a positive solution to (6-2). Then for any p > 1, 

(6.7) (p - 1) ^_ (! ,_ 1)? <i. 

The result also holds on noncompact manifolds assuming the the left hand side 
is bounded. Note that when p — 2, the above reduces to Li-Yau's estimate. Hence 
the corollary provides another nonlinear analogue of Li-Yau's estimate for the heat 
equation. 
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Remark 6.4. It is easy to check that for p > 1, the function 

. v I -nil, .p-i .s-n/ P T(n/2 + 1) / |x-a;o| 

(6.8) w = 7T ' (p * pt) /p > - ■ exp ' 



r(n/p* + l) ^ (ip*p-ip)iAr 

is a fundamental solution of (6.2) on R n , and achieves equality in the estimate (6.7). 
This demonstrates the sharpness of the estimate. Here p* = -^-^ . 

Remark 6.5. It seems that an estimate in the above sharp form is not known 
except in dimension one [EV1]. It should be useful in the study of the regularity of 
weak solutions of (6.2). 

When p = 1 we have the following result. 

Theorem 6.6. Let M be a compact manifold with nonnegative Ricci curvature. 
Let u be a smooth solution to (6.1). Then 

(6.9) \^ u \ +Ut <!Lzl. 

Since the equation becomes degenerate when Vn = 0, the above computation 
needs extra justification. First we introduce the notion of a weak solution to as in 
[He]. For any compact subset K, we define the functional 

J p (w;K) = [ -\\7w\ p + w{\Vu\ p + ut)dfi 

JKx{t} P 

for w G C = {ip G C°(M x [0,oo)) | p(-,t) G C" M (M), and for every x G M,ip(x,-) G 
C ; l ^(0, oo)} such that {w ^ u} CC M. A function u is called a continuous weak 
solution to (6.1) if 

Jl(u;K)<JZ(w;K) 

for any w G C. Note that when u is smooth (locally), satisfying (6.1) with |Vu| / 0, 
we have that 

Jl{u;K)-Jl{w;K) = f -(\X7u\ p - \X7w\ p ) + (u - w) divfl V«| p - 2 Vu) dix 

Jk P 

< [ -(\Vu\ p -\Vw\ p )-\X7u\ p + \X7u\ p - 1 \X7w\dn 
Jk P 

< 

by Young's inequality. To obtain weak solutions more generally, we use an e-regularization 
process, replacing the equation (6.1) by the approximate version 

(6.10) ^ = A p , e (ue) = div(/f/ 2 - 1 V Ue ) - f p/2 

where f e = |Vw e | +e. Since (6.10) is strictly parabolic, one can apply the established 
theory of parabolic equations to obtain a solution u e and take the limit e — > to 
obtain a continuous weak solution to (6.1) as in [He]. We next show that an estimate 
analogous to (in fact, somewhat stronger than) (6.9) can be obtained for solutions 
u € to (6.10), and with this estimate recover (6.9) in the limit as u € — * u. 
It will be convenient to introduce the notation 

C^) = DA p , E [u E ](V) = div(/ e p / 2 -U e W) -p/f /a_1 <Vu e , VV) 
for the linearization of the operator A p>£ , where 

^=id+(p-2) Vt -f Vl - . 

It is easy to check that 

dt dt ' 
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and a computation as before establishes the following Bochner-type formula. 
Lemma 6.7. 

Ce{fe) = 2f v J 2 - x (|VVu E | 2 + Ric(Vu E , Vlie)) + 2<Vu E , V(A p , e (u e ))> 
(6.11) +(|-l)!V/ E | 2 /f /2 " 2 - 



In particular, 



(c t - f? /2 = P fr 2 (Ric(V We , Vu E ) + |VVu«ft. 



Thus we obtain the analog of Corollary 6.2 for solutions to our approximate equa- 
tion (6.10). 

Lemma 6.8. For all a > 0, 

satisfies 

(6.12) (c e - F a , e = p/r 2 (Ric(Vu E , Vu e ) + | VVu e \ 2 Af ) . 

In the same way as before, we note 

Fi, e = div(/ E p/2 - 1 Vn E ) = fj 2 - 1 %r At (VV Ue ) 



so 



/f- 2 |VV«H > iF 2 E , 

which, with the above lemma, implies the following global estimate. 

Theorem 6.9. Suppose u t is a positive solution to (6.10) on the compact manifold 
M of non-negative Ricci curvature. Then, for all p > 1 and all t > 0, 

(6 . 13) ( |v„,|* +e) - + ^<i 

Hence Theorem 6.3 and Theorem 6.6 hold for the weak solutions obtained via the 
e-regularization process. 

The Bochner-type formula (6.6) is effective enough to give a nonlinear entropy 
formula (see [N] for the entropy formula for the heat equation). For this purpose, 
we first observe that f u p_1 dfx is preserved under the equation (6.2). A little less 
obvious, perhaps, is the following conservation law. 

Proposition 6.10. For any smooth yj, 

( 6.14) ± lj^^l u ((»-/:)*).-'* 

In particular, if ip satisfies 

(6.15) (1-A^ = 



then f M ipv p ~ dfj, is a constant along the equation (6-2) and (6.15). 
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Proof. Notice that Vv = -^Vm, hence (p-l)^ 1 1 Vv\ p - 2 Vv = -v v - 1 \ Vw| p " 2 Vw. 
Direct computation shows that 

f [ ^ p - 1 d / u = / ^r'i^/- 1 ), 

"* J M J M 

-(p-1) / |Vw| p_2 (Vw,VV')v p - 1 . 
Jm 

The result follows from the observations that Vv p_1 = — u p_1 Vif and 
-/ \\7u\ p - 2 (A(\7i>),\7v p - 1 ) = -[ \Vu\ p - 2 (Vip,Vv p - 1 ) 

J M J M 

-(p-2) / |V«| p " 4 (VV',Vk}(Vk,V«^ 1 ) 
= (p-1) / |Vu| p -V _1 (Vi/>,Vu). 

JM 

q.e.d. 

Remark 6.11. In the above we do not really need that ip is smooth. The argument 
carries over assuming only the summability of the integrands. 

Let (dij ) be the inverse of (A lj ) . Explicitly, = g%j — P ^f — j 2, > an d can be viewed 
as a metric tensor. Let v be a positive solution to (6.2) satisfying J M v p ~ 1 dfi = 1. 
Define 

J\f p (v,t) = / v p ~ 1 ud{i ~ — logi 
Jm P 

and denote the first term on the right hand side by N(v, t). Note that u = — log(t; p_ ). 
Following [N], (see also [FIN]), define F p (v, t) = f t M p {v,t) and W p (v,t) = f t (tN p {v,t)). 
More explicitly, motivated by (6.8), if we write 

v*- 1 = 1 ^ r(n/2 + 1) e-J ^ 

7r n/2( p »p-l p )f T(n/p*+l) ' 

then define 

W p {<p,t)= f {t\X7<p\ p + <p - n^"- 1 dp. 
Jm 

For the case p = 1, it is helpful to write W in terms of u, 

W P (u,t) = f (t\Vu\ p + u--logt + \og( w 3 n {\ +l \ \ -n) e~ u dix, 

which becomes 

Wi (u, t) = I (t\Vu\ + u - n log t + log V(jl ^ - - n) e" 11 d M 

when p = 1. Note that when p = 2, this is precisely the entropy defined in [N]. 
For the entropy quantity W p (v, t) we have the following nonlinear entropy formula. 

Theorem 6.12. Let v be a positive solution to (6.2) satisfying f M v p ~ 1 dfi — 1. 
Then 

(6.16) ^-W p (v,t) = -tp f ( \f p/2 - 1 V l V j u- — aij + f p - 2 R lJ u i uAv p - 1 dp,. 
dt Jm \\ f P a 

When p — 2, this recovers the entropy formula of [N]. 
For the proof we need the following result. 
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Proposition 6.13. 

(6.17) 4 MM) = / F 1 v p - 1 d f i= [ f' 2 ' 1 trACVVuK" 1 dp, 
dt J M J M 

(6.18) / Fiv^dn = -p f f p ~ 2 (\VVu\ 2 A +R l ,u t u J )v p - 1 d f i. 
dt J M J M 

Proof, (of Proposition.) Direct calculation shows that 

^ - £^ w = (p - l)\Vu\ p - (p - 2) div(|Vu| p - 2 Vu). 

Hence the first identity of Proposition 6.13 follows from (6.14) together with the 
observation 

f F 1 v p - 1 = f div(\X7u\ p - 2 X7u)v p - 1 = f |Vu|V _1 . 

JM JM JM 

The second identity of Proposition 6.13 follows from (6.6) and (6.14). q.e.d. 

Note that in the case that the Ricci curvature of M is nonnegative, (6.18) yields the 
monotonicity of 'energy' 



F{v,t)= [ |Vu|V -1 . 

JM 



Now Theorem 6.12 follows from the above proposition similarly as in [N] (also [FIN]), 
by observing that = t-^f- + 2T V and completing the square. This point of view 
is taken from physics [Ev, P]. The case p = 1 (Theorem 1.4) can be shown similarly. 

By ODE considerations and the Cauchy-Schwarz inequality, we also have the 
following result. 

Corollary 6.14. Assume that M has nonnegative Ricci curvature. Then 

In particular, any positive ancient solution to (6.2) must be a constant. 

Summarizing we have that if M has nonnegative Ricci curvature then J\f p (v,t) is 
a monotone non-increasing concave function in logt. 

To deal with the potential vanishing of Vu, we work via the approximation scheme 
as before. First, (6.17) can be justified in view of Remark 6.11. For (6.18), we need 
to appeal to an approximation argument. Let w e be a solution to (6.10) and define 
tfit similarly as before. We have the following pointwise computations. 

Proposition 6.15. Let u e , tp E and / E be functions defined as before. Let 

V E = 2div(|Vu E | p " 2 Vit E ) - |Vu E | p , W e = t(2div(|Vv9 e | p ~ 2 V<y2 E ) - \\7ip € \ p ) + ip € - n. 

Then 

(6.20)^ - £?j V e = -pfr 2 (|VVu e |i + Ric(Vu E , V« e )) 



+ Ric(Vw E , Vu £ ) 



tp 

(6.21) +{p - 2) (| Vit E | p - div(|Vu E | p ~ 2 Vu E )) . 

Notice that 

/ V ( v p - 1 =2[ rVu e ,Vu)\Vu e \ p - 1 v p - 1 - IVw^V 1 ^ / \Vu\ p v p - 1 =T 

JM JM JM 



GRADIENT ESTIMATES AND AN ENTROPY FORMULA 



27, 



and 



4 / = ~P [ /r 2 (|VVu E |^+Ric(Vu e ,V« e ))« P 

at J M J M 



J 



/f/ 2 - 1 f <w e , w 1 ) + ( P -2) (v ^; VK) (Vn £ ,v^- 1 ) 

-P / /f /2 " 1 (V«e,VK)« p - 1 +f /"^-V-^Vu.VV.) 
->• -p / f~ 2 (|VVw|i + Ric(Vu, Vu))/" 1 

JM 

as e — > 0. Hence (6.18) can be established for weak solutions which can be properly 
approximated by the regularization process. 

Strictly speaking, the above prool can be applied only for the case when M is 
compact. When M is non-compact, lurther work is needed to justily the integrations 
by parts involved. On complete noncompact manifolds with nonnegative Ricci cur- 
vature (or for which the Ricci curvature is bounded from below) the entropy formula 
has been established rigorously in the forthcoming book [Chow, et al p2] for the 
case p = 2, with full justification for such manipulations. The argument there can 
be adapted to the p > 1 case. Since this is rather technical we refer the readers to 
[Chow, et al p2] for details. 

In [N], the entropy formula for the heat equation was used to show that a man- 
ifold with nonnegative Ricci curvature and sharp logarithmic Sobolev inequality 
must be isometric to R n . It is interesting to study the relation between the L p - 
logarithmic Sobolev inequality and the entropy formula (6.16). First recall the sharp 
L p -logarithmic Sobolev inequality (see for example [DpDG]). 

Theorem 6.16. Let M be Euclidean space R n . Then for any p > 1, and any w 
with J \w\ p dx — 1, 

(6.22) j \w\ v log \w\ v dx<^ log (c p , n j \Vw\ p dx^j 



where 



^ n =£r^r 1 ^^ r(n/2+i)N " " 



n\ej \T{n/p* + 1) 



The above inequality with sharp constant was established after the work of [DD, 
G] (see also [E]), and has the following form connected with the entropy quantity. 
(For the L -version, see [Bo] as well as [Be].) 

Proposition 6.17. If on (M,g), (6.22) holds then 

W p {tp,t) > 

for 

p-i _ 1 T(n/2 + l) e"" 



7r«/2(p*p-ip)f F(n/p* + 1) fp 
with J M v p ~ 1 = 1, where t is just a scaling factor. 

Proof. The proof is a technically more complicated version of the p = 2 case 
shown in [Chow, et al pi], pages 247-249. We leave the detailed checking to the 
interested reader. q.e.d. 

We shall next prove the following characterization of R™ among manifolds with 
nonnegative Ricci curvature by the sharp L p -logarithmic Sobolev inequality. 
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Theorem 6.18. Let (M,g) be a complete, connected Riemannian manifold with 
nonnegative Ricci curvature. Assume that (6.22) holds for some p > 1 with the sharp 
Then M is isometric to R n . 

Proof. By an argument along the lines of [P], Section 3, one can show that for the 
fundamental solution H to (6.2), linit^o W P (H, t) — 0. Then, the entropy formula 
ensures that W P (H, t) < 0. However, in the presence of a sharp L p -logarithmic 
Sobolev inequality, one may prove, as in Corollary 6.17, that W P (H, t) > for t > 0. 
Hence we can conclude that -^W P (H,t) = 0. Writing u = — (p — l)\ogH, and 
/ = |Vtt| 2 we have that 

(6 ' 23) ^ = ^i s ^FTH 2 

off of the set of critical points of u, which we shall denote by E. 

Equation (6.23) implies in particular that VVu > on M \ E, so VVw > on all 
of M. It follows, then, that p G E if and only if u(p) = inf 96 M u(q) = in. We shall 
prove later that E consists of precisely one point; for the present we note that these 
facts imply at least that E is a connected (in fact, convex) subset of M. Indeed, if 
Pi, P2 £ E, and 7 : [0, 1] — > M™ is a unit-speed geodesic with 7(0) = pi and 7(Z) = P2, 
then we have 

d 2 

— u( 7 (a)) = (V^u) ( 7 (a)). 

Since ^(7(5)) > u(-y(0)) = u(-y(l)) = in, the weak convexity of u implies u(-y(s)) = rn 
and so j([0,l]) C E. 

Let us now consider the structure of an arbitrary connected component Mi of 
M \ E. We may write Mi w I x N for an interval J and a hypersurface N defined by 

N = Mi n {u = u(p )} 

for some arbitrary po € Mi. Now we choose a orthonormal frame such that ei = 
Then the equality (6.23) implies that 



Mil 



tp(p- l)u\ 2 ' 

«l/3 = 0, 

1 



U. 



a: a 



tpu\ 

for 2 < q,/3 < n. This shows that |Vw| is constant along the level set of u. We 
shall use the above information to find the explicit form of the metric. Recall that 
u a p = h a @Ui, where h a p is the second fundamental form of the level set hypersurface 
of u. Let U — p|Vm| p_1 . It is easy to see that the mean curvature of this hypersurface 



is H= 2=±. Hence 



n — 1 

Ui = 1, AU = Un + HUi = 



u ■ 

If we fix a level set and parametrize the other ones using s, the oriented distance 
from this fixed one, we have that 

1 



p(p - 1)U% 

which implies that (u P_1 ) s = i. Thus 



(6.2.1) in = [ ^ +B 
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where B is the value of |Vu| p 1 on the originally fixed hypersurface. This implies 
that the level hypersurfaces at distance s have the second fundamental form 

, 1 

This further implies that the metric on N x / is of the form 

(6.25) g i3 = ds 2 + (s + pB) 2 g N 

where qn is the metric of the hypersurface N. Since M is complete, and Mi a 
connected component of M \ E, we see that I = (— pB, oo). 

Now, from equation (6.25), one may compute that Km{X, Vm, X, V«) = along 
N for any X tangential to N and hence (either again from (6.25) or from the Gauss 
equations) that 

(6.26) Rc N (X,X) = Rc M (X,X) + (n - l)g N (X,X). 

Since the Ricci curvature of M is assumed to be non-negative, it follows immediately 
from Myer's theorem that N is compact and, moreover, that the diameter of the 
hypersurface N s of signed distance s £ I from N satisfies the bound 

(6.27) diam(jV s ) < rr(s + pB). 

Thus lim s __ P B diam(iV s ) = and dM\ must consist of a single point O £ E. In 
fact, we must have M = Mi U {0} as otherwise the removal of O would disconnect 
M - an impossibility if dim(M) > 2. Thus M = Mi U {O} in this case. (Note that 
since the foregoing considerations imply that M must be non-compact, the theorem 
has already been established in the case dim(M) = 1.) So N (and hence all N s ) 
are geodesic spheres about the point O, and the representation (6.25) on M \ {O} 
implies that M is isometric to R". q.e.d. 

Remark 6.19. Here we assumed the existence and regularity of the fundamental 
solution, together with the finiteness of all the integrals involved in the entropy. 
These assumptions certainly require justification. 



7. Localization 

In this section, carrying on in the notation of the last, we localize the estimate (6.7) 
of Theorem 6.3. Here, in considering the parabolic equation (6.1) we encounter some 
technical complications beyond those encountered in Section 2 for the derivation of 
the local estimate for the elliptic version, including additional complications caused 
by the nonlinearity /degeneracy of the equation. 

First we focus on the case 1 < p < 2. As in Section 2, assume that u is a solution 
on B(xo, R) and let 8{t) : [0, oo) — » [0, 1] be a cut-off function satisfying < 10 

and 0" > -100 > -10. Now let rj(x) = P (^). It is easy to see that 

Tj P 

The proof of Lemma 2.2 implies the following estimate. 

Lemma 7.1. Assume that the sectional curvature is bounded from, below on 
B(x Q ,R) by -K 2 . Then 

(7.2) At? +(p-2)I^>- Ci y fl) 
where d{n,KR) = 40(n +p - 2)(1 + KB) - 20. 
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Theorem 7.2. Assume that the sectional curvature of M on B(xo,R) satisfies 
Km > —K 2 and u is a solution of (6.1 ). Assume further that 1 < p < 2 and u t < 0. 
Then for any a > 1, on B(xo, -§)> 
(7.3) 

(%na 2 (^ + ^y , (4n(n - l)K 2 ) P ' 2 

Remark 7.3. The above result is a parabolic analogue of Theorem 1.1. Indeed, 
if it is a solution to the elliptic problem, we have ut = 0, hence one can apply the 
theorem and recover the elliptic estimate (albeit with a worse constant). 

Proof. Clearly, tF a n achieves a maximum over B(xo,R) x [0, T] somewhere. For 
the purpose of bounding F a from above, we may assume that r\F a > at this max- 
imum point (j/o,io); this implies in particular that io > 0. Applying the maximum 
principle at (j/o, to), we have that ( Jj — C) (tF a n) > 0. 

Direct computation shows that 

(jt-cj(tF a ri) = toi(jt-c) F a +tF a (^- t -C S jn-2tf ,2 - 1 (Vn,VF a ) A 
+riF a . 

Here (•, -)a denotes the product with respect to A lj . Using that V(rjF a ) = at the 
maximum point, with formula (6.6) and the simple estimate on V ^ A from Section 
2 we find that 

< (| -*)(**,,,) 

< -Ptvf- 2 (|VV«ft - (n - l)K 2 f) + tF a (jft-Cjr, 

+ 2tr/^F a ^+ v F a . 
V 

We first derive an effective estimate on F a (Jj — C) r\. Direct calculation using (6.1) 
yields 

- (a v + {p- 2)Mi^ fl^F a +pf ,2 - 1 F a (Vu, V n ) 
+{ 2-p)F^^L. 
We shall estimate each term in turn. Define A %i = g lj — ^y- . It is easy to see that 

< (2-p)F a f^ 2 - 2 (\7 V ,\7f) A 

< (2-p)F a ft 2 - 2 \Vn\\Vf\ A . 

Choose a normal frame so that e\ — • Then Uj = for all j > 2, and f k — 2u\ h u\ 
for k > 1. Hence 

F a .r /2 - 2 wf\ A = 2 Fa .r' 2 - 2 = 2^/ p/2 - 3/2 ■ 
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On the other hand, 

IVVuji = u lk + 2{p - 2)4i + (p - 2) 2 Wll 

= (p - l) 2 f 11 + 2(p - 1) X] "i 1 + u i fe 

j>2 j,fc>2 
2 



> 2(p - 1) J] «?i + - (tTA(VVu)) 



J>2 



Hence we have that 



2 p — 1 f] j In 

Lemma 7.1 implies that 

(7.5) II = - (a„ + (p- 2)»^) Z'/ 2 - 1 ^ < CiW) r / 2 -i Fa . 
Furthermore, it is easy to see that 

(7.6) /// = p/ p/2 - 1 F a (V«,V??) ^p/^i^V??! 
and 

(7.7) = (2-p)F 1 F a (VM / V?7) < (2-p)FiF a |V77|r 1/2 . 

Since / p_2 |VVu|^ > we split ??F a = y - az with y = r)\Vu\ v and 2 = rju t . 

Also we can write rfF? as 



(7.8) — :// azf | 2—^1;, | I - -') y 



Note y — z — — az) + (1 — ^)y > 0. Combining (7.4)— (7.7), we arrive at 
< -^^/^IVVnlU^-W/^^ + ^IV^F 2 /- 1 
+ C 1 (n,gfl) 7?r/2 _ 1 ^ +y/£? i Fa|v??|?? + (2 _ p)FlFa | Vl?h/ -i/2 

fP/2-l/r_|T7„|2., J ? 2 - F « 



+/ p/ ^ 1 F a |Vr;r + 

3p-2 2 2 _a 2 2-p 2 _a IV^I 2 

< 7T— (y-z) +p(n-l)y pK +——(y-az)y p J— V 

ZT1 P 1 Tj P 

CAn,KR) , , !_p , , i_i |Vr?| 

It fj p 

+(2 - p)(y - ,)(y - a*)*"* + (y - a V* ^ + 

jj p 7] p 

1 2 a — 1 /a — l\ 2 

+2— (_-) y 

+2(n - l)y 2 -t^ 2 + (y - a*) V* + ^^(y - -)/"§ 

20 !_i (2-p) 2 _i (2-pXa-l) !_. 
+ -^-(y-a«)y p + — (y - az ) V ^ (y-az)y 

100. . i_i y — az 

+^-(y~ az )y p + — i — • 
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The above gives a local estimate if we assume u t < (namely z > 0). (Assuming 
«t < initially, this can be established by the maximum principle, in view of (6.4).) 
In this case, we have that y — az < y. 

We make use of the above computation as follows. If we do not have 

/ \ 2 

(7.9) y 2/p < 4n(n - 1)K 2 



which implies that y — az is bounded by the same bound, we can drop the terms only 
involving the powers of y to still have the last estimate on the functions of y — az. 
By the elementary inequality 2ab < a 2 + b 2 , we have that 



- , , _ Ci(n,KR) + 100 r_ 

(y-az) [ -—z-y + ^ y " + -j^ y 



\ na- 1 

I 90 4- O. - n\ — 

a J 4(a - 1)1 R 2 
< ±(c 1 + 100 + (20 + (2 - p) £^1) 2 ^ ) (, - a,)- f . 



< (y - az) (d + 100 + (20 + (2 - p)^) ' 77 ^- TT ) -^y 1 



Here we have used y — az < y. This takes care of the terms in which (y — az) appears 
linearly. 

Making use of y — az < y again, we have that 
< Ay -az) +-^{y-az) ? H —(y-az) p + 



2na 2Ka ' R 2 a ' aR Ka ' t 

with 

C 2 = Ci(n, tffl) + 100 + (20 + (2 - p )^±) 2 na \ + 50^. 

a i(a — 1) p — 1 

Another use of the above elementary inequality reduces it to 

(7.10) < - JL(„ - a,) 2 + g(y - a,) 2 -§ + ^ 

with C2 = C 2 + n(2 - p) 2 . Using Young's inequality a& < ^- + with ^ + = 1, 
we finally have that 

(7.11) 0<-^ { y-azf + (% + ^) { y-az?-i 
where 



C 3 = f (4na 2 (2-p))" 
From (7.11) it is then easy to derive the theorem. q.e.d. 
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